Phase Transition from QMC Hyperonic Matter to Deconfined Quark Matter 



J. D. Carroll0 D. B. Lcinwcbcr. and A. G. Williams 

Centre for the Subatomic Structure of Matter (CSSM), 
Department of Physics, University of Adelaide, SA 5005, Australia 



o 
o 

(N 
u , 

< 



> 
00 



A. W. Thomas 

Thomas Jefferson National Accelerator Facility, 12000 Jefferson Ave., Newport News, VA 23606, USA 
College of William and Mary, Williamsburg, VA 23187, USA. and 

Centre for the Subatomic Structure of Matter (CSSM), 
Department of Physics, University of Adelaide, SA 5005, Australia 
(Dated: April 29, 2009) 

We investigate the possibility and consequences of phase transitions from an equation of state 
(EoS) describing nucleons and hyperons interacting via mean-fields of a, uj, and p mesons in the 
recently improved Quark-Meson Coupling (QMC) model to an EoS describing a Fermi gas of quarks 
in an MIT bag. The transition to a mixed phase of baryons and deconfined quarks, and subsequently 
to a pure deconfined quark phase is described using the method of Glendenning. The overall EoS for 
the three phases is calculated for various scenarios and these are used to calculate stellar solutions 
using the Tolman-Oppenheimer-Volkoff equations. The results are compared to recent experimental 
data and the validity of each case is discussed with consequences for determining the species content 
of the interior of neutron stars. 

PACS numbers: 26.60.Kp, 21.65.Qr, 12.39.-x 



O ■ I. INTRODUCTION 



The use of liadronic models to describe high density matter enables us to investigate both the microscopic world of 
atomic nuclei and the macroscopic world of compact stellar objects, encompassing an enormous range of scales. The 
results of these investigations provide deep fundamental insight into the way that the world around us is constructed. 
Experimental data from both extremes of scale aid in constraining such models, from the saturation properties 
\^ I of nuclear matter to the observed properties of neutron stars [H, 0, Q . The literature 0, H, 0, 0, H, provides a 
T— I ■ plethora of models for the EoS of hadronic matter, at least some of which have been successfully applied to calculate 

■ the properties of finite nuclei. There are also important constraints from data involving heavy-ion collisions [lol. [ll]|. 
. Many of these EoS have also been applied to neutron star features as well. However, the amount of data available for 

■ neutron stars (or com pac t stellar objects) is very limited, with only a single result containing both a mass and radius 
OO result simultaneously |12| | and even that has been recently disputed [l^l ■ 

■ With such a lack of constraining data, our focus shifts to finding models which better reflect the physics that is 
^ I expected to be important under the conditions which we are investigating. A prime example of this is that at the 

densities which we consider to be interesting for this investigation (1-10 times nuclear density, po ~ 0.16 fm^""*) it is 
possible that either hyperonic matter (in which strangeness carrying baryons become energetically favourable as the 
Fermi sea fills), quark matter (in which it becomes energetically favourable that the quarks inside the baryons become 
deconfined) or a mixed phase of these is present, rather than the more traditional treatment of nucleons alone. 

We construct a model of high density matter which is globally charge neutral, color neutral, rotationally symmetric, 
and in a state that is energetically favourable. For this purpose we consider hadronic matter modelled by the Quark- 
Meson Coupling (QMC) model [3, Ell, which was recently improved through the self-consistent inclusion of color 
hyperfine interactions [l6[. While this improvement had no significant effect on the binding of nucleons, it led to 
impressive results for finite hypernuclei [l3j- We follow the method of Glendenning [3 to produce a mixed phase of 
hyperonic matter and deconfined quark matter under total mechanical stability, then a pure deconfined quark matter 
phase with relativistic non- interacting quarks. 

We begin with a brief presentation of Relativistic Mean-Field Theory in Section [IT] to establish a foundation with 
which to discuss the general formalism for the QMC model, including new additions, in Section Hill Deconfined quark 
matter is discussed in Section IIVI This is followed by Section |V] providing a summary of the requirements for and 
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method to construct a phase transition from a hadronic phase to a mixed phase and from a mixed phase to a quark 
phase. Stellar solutions are calculated in Section IVII and a summary of our results is presented in Section IVIII with 
conclusions in Section IVIIII 



II. RELATIVISTIC MEAN-FIELD THEORY 



We introduce the mean-field description of nuclear matter using the classic example of Quantum Hadrodynamics 
(QHD) [l^, [23|. Although the Quark-Meson Coupling model (QMC) has a fundamentally different starting point, 
namely the self-consistent modification of the structure of a hadron immersed in the nuclear medium [2ll . [U[23| , in 
practice the equations for nuclear matter involve only a few key differences. We summarize those in the next section. 
The original formulation of QHD included only nucleons interacting with scalar-isoscalar, cr, and vector-isoscalar, uj, 
mesons. This was later expanded to include the vector-isovector, p, and subsequently the entire octet of baryons, 
B G {p, n, A, S+, S~, S°, 2~}, with global charge neutrality upheld via leptons, £ S {e~ 

The Lagrangian density for QHD is 



C = 



k 



i'k - gukuj'' - gpf(k) ■ P^) - [Mk - gak(j)\ ipk 



(1) 



where the index k G {iV, A, S, 
spinors for these 



2 F '2 

1} represents each isospin group of the baryon states, and tpk corresponds to the Dirac 



N 




The vector field tensors are 

The third components of the isospin matrices are 



T{N)3 



1 

-1 



r(A)3 



[0 ], 



'-(S)3 






-1 



(2) 



(3) 



(4) 



^pl is a spinor for the lepton states, and SC are renormalisation terms. We do not include pions here as they provide 
no contribution to the mean-field, because the ground state of nuclear matter is parity-even. We have neglected 
nonlinear meson terms in this description for comparison purposes, though it has been shown that the inclusion of 
nonlinear scalar meson terms produces a framework consistent with the QMC model without the added hyperfine 
interaction [2^ . The values of the baryon and meson masses in vacuum are summarized in Table HI 

Assuming that the baryon density is sufficiently large, we use a Mean-Field Approximation (MFA) with physical 
parameters (breaking charge symmetry) in which the meson fields are replaced by their classical vacuum expectation 
values. With this condition, the renormalisation terms can be neglected. 

By enforcing rotational symmetry and working in the frame where the matter as a whole is at rest, we set all of 
the 3-vector components of the vector meson fields to zero, leaving only the temporal components. Furthermore, by 



TABLE L The vacuum (physical) baryon and meson masses (in units of MeV) as used here [25l |. 





M„ 


Ma 




Mj20 




Ms- 


Mho 


938.27 


939.57 


1115.68 


1197.45 


1192.64 


1189.37 


1321.31 


1314.83 














mp 






550.0 






782.6 




775.8 
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enforcing isospin symmetry we remove all charged meson states. Consequently, because the mean-fields are constant, 
all meson derivative terms vanish, and thus so do the vector field tensors. The only non-zero components of the 
vector meson mean fields are then the time components, (w^) ~ {uj)S^'^ and (p^) — (p)(5'^°. Similarly, only the third 
component of the p meson mean field in iso-space is non-zero, corresponding to the uncharged p meson. 

The couplings of the mesons to the baryons are found via SU(6) flavor-symmetry [2^. This produces the following 
relations for the a and lo couplings to each isospin group (and hence each baryon B in that isospin group) 

111 111 

Using the formalism as above with isospin expressed explicitly in the Lagrangian density, the couplings of the p meson 
to the octet baryons are unified, thus by specifying g^N ^ g^N, and gp we are therefore able to determine the couplings 
to the remaining baryons. 

By evaluating the equations of motion from the Euler-Lagrangc equations 

dC dC 

- Op = 0, (6) 

we find the mean-field equations for each of the mesons, as well as the baryons. The equations for the meson fields 
are 

^ — ' ml 

B 

(-) = E ^(^b7Vb) = E '-^(^Ub), (8) 
^ — ' m^, ^ — ' mf, 

g ^ B ^ 

(P^ - E ^(V5.7"r(.)3^.) - E Zi^i^hm^^) - E 7^«I^b^b), (9) 

k P k P BP 

where the sum over B corresponds to the sum over the octet baryon states, and the sum over k corresponds to the 
sum over isospin groups. I3B is the third component of isospin of baryon B, as found in the diagonal elements of T(^k)3- 
(lu), (p), and (a) are proportional to the conserved baryon density, isospin density and scalar density respectively, 
where the scalar density is calculated self-consistently. 

The Euler-Lagrangc equations also provide a Dirac equation for the baryons 

E ^ - .9c.B7°(^) - gpl°hB{p) - Mb + gaBifj)] ^B = 0. (10) 

B 

At this point, we can define the baryon effective mass as 

M*B^MB-gaB{o), (11) 

and the baryon chemical potential (also known as the Fermi energy, the energy associated with the Dirac equation) 
as 



^kl^ + [Ar^f + g^u) + gphBip). (12) 
The chemical potentials for the leptons are found via 



pi = ^k%+ml. (13) 

The energy density, f , and pressure, P, for the EoS can be obtained using the relations for the energy-momentum 
tensor (where is the 4- velocity) 

(T^'') = {£ + P) u^'u'' + Pg^^ => P= -{T"), £ ^ (r"°), (14) 

3 

since = and uqvP — —1, where g^'^ here is the inverse metric tensor having a negative temporal component, 
and T^" is the energy-momentum tensor. In accordance with Noether's Theorem, the relation between the energy 
momentum tensor and the Lagrangian density is 

T-.-g-/: + a^^^, (15) 
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and we find the Hartree-level energy density and pressure for the system as a sum of contributions from baryons, B; 
leptons, £; and mesons, m to be 

j=B,l.m 



^ = E 

j=B.f.m 



(16) 



v-^ (2J, + 1) f e{kF, - \k\) V- 1 2 / \2 1 2/ \2 

3(2^)3 J ^W^iMfy ^i^J ^ 2 ^ ' 

(17) 

where Ji is the spin of particle i (J,; = ^ V i G {B,£}) which in this case accounts for the availability of both up 
and down spin-states. 9{x) is the Hcaviside Step Function. Note that the pressure arising from the vector mesons is 
positive, while it is negative for the scalar meson. 
The total baryon density, p, can be calculated via 

(2Ji3 + 1) 

B B 



- E/^^ = E (2^)3 j (^(kFs - |fc|) d% (18) 



where in symmetric matter, the Fermi momenta are related via kp = kp^ = kpp, and the binding energy per baryon, 
E. is determined via 



E : 



l(^£-Y,MsPB 



(19) 



The couplings g^N and g^jN are determined such that symmetric nuclear matter (in which Pp = Pn = 0.5p) satu- 
rates with the appropriate minimum in the binding energy per baryon of Eq = —15.86 MeV at a nuclear density of 
po = 0.16 fm~'^. The couplings for QHD which provide a fit to saturated nuclear matter are shown in Table HIl 
The coupling gp is fixed such that the nucleon symmetry energy, given by 

1 fr2 1 1,2 



12.2m^ ^ - 12 ^kl + iM;r 12 VWTim' 

is reproduced at saturation as (asyni)o = 32.5 MeV. 

The chemical potential for any particle, pi, can be related to two independent chemical potentials — we choose 
that of the neutron {pn) and the electron (pe) — and thus we use a general relation 

p^ ^ B,pn - Q^pe ; i £ {p, n. A, E+, s", S" , S°, S" , ^}, (21) 

where Bi and Qi are the baryon (unitless) and electric (in units of the proton charge) charges respectively. For 
example, the proton has Bp = -1-1 and Qp = +1, so it must satisfy Pp = Pn — A*e which is familiar as /3-equilibrium. 
Since neutrinos are able to escape the star, we consider p^ = 0. Leptons have Bi = 0, and all baryons have Bb — +1. 
The relations between the chemical potentials are therefore derived to be 

MA = Ms" = Ms" = Mm 

Ms- = Ms- = A'n + Me: ^22) 
Mp ^ MS+ = Mn ^ Me I 

Mm = Me- 

The EoS for QHD can be obtained by finding solutions to Eqs. ([7H9]) subject to charge neutrality, conservation of a 
chosen total baryon number, and equivalence of chemical potentials. These conditions can be summarised as 

P ^E^B^p^ \ i G {p, n. A, S+ , , , S° , , £} . (23) 
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With these conditions, we are able to find the EoS for QHD. It should be noted that, as with many relativistic models 
for baryonic matter, once we include more than one species of baryon this model eventually produces baryons with 
negative effective masses at sufficiently high densities (p > 1 fm~^). This is a direct result of the linear nature of the 
effective mass as shown in Eq. (|lip . As the Fermi energy (sec Eq. (|12p ) approaches zero, the cost associated with 
producing baryon-anti-baryon pairs is reduced and at this point the model breaks down. From a more physical point 
of view, as the density rises one would expect that the internal structure of the baryons should play a role in the 
dynamics. Indeed, within the QMC model, the response of the internal structure of the baryons to the applied mean 
scalar field ensures that no baryon mass ever becomes negative. We now describe the essential changes associated 
with the QMC model. 



TABLE II: Couplings for QHD with the octet of baryons, fit to saturation of nuclear matter. 







9p 


10.644 


13.179 


6.976 



III. QMC MODEL 

Like QHD, QMC is a relativistic quantum field theory formulated in terms of the exchange of scalar and vector 
mesons. However, in contrast with QHD these mesons couple not to structureless baryons but to clusters of confined 
quarks. As the density of the medium grows and the mean scalar and vector fields grow, the structure of the clusters 
adjusts self-consistently in response to the mean-field coupling. While such a model would be extremely complicated to 
solve in general, it has been shown by Guichon et al. [27t that in finite nuclei one should expect the Born-Oppenheimer 
approximation to be good at the 3% level. Of course, in nuclear matter it is exact at mean-field level. 

Within the Born-Oppenheimer approximation, the major effect of including the structure of the baryon is that the 
internal quark wave functions respond in a way that opposes the applied scalar field. To a very good approximation 
this physics is described through the "scalar polarizability," d, which in analo gy w ith the electric polarizability 
describes the term in the baryon effective mass quadratic in the apphed scalar field [2 ll. [28l. [29l. [sol. [sH . Recent explicit 
calculations of the equivalent energy functional for the QMC model have demonstrated the very natural link between 
the existence of the scalar polarizability and the many-body forces, or equivalently the density dependence, associated 
with successful, phenomcnological forces of the Skyrmc type [T^ . [Tsj . In nuclear matter the scalar polarizability is 
the only effect of the internal structure in mean-field approximation. On the other hand, in finite nuclei the variation 
of the vector field across the hadronic volume also leads to a spin-orbit term in the nucleon energy (27j . 

Once one chooses a quark model for the baryons, and specifics the quark-level meson couplings, there are no new 
parameters associated with introducing any species of baryon into the nuclear matter. Given the well known lack of 
experimental constraints on the forces between nucleons and hyperons, let alone hyperons and hyperons, which will 
be of great practical importance as the nuclear density rises above (2~3)/7oj this is a particularly attractive feature 
of the QMC approach and it is crucial for our current investigation. Indeed, we point to the very exciting recent 
results of the QMC model, modified to include the effect of the scalar field on the hyperfine interaction, which led to 
A hypernuclei being bound in quite good agreement with experiment and E hypernuclei being unbound because of 
the modification of the hyperfine interaction [l7t - thus yielding a very natural explanation of this observed fact. We 
note the success that this description has found for finite nuclei as noted in [Tsj . 

While we focus on the MIT bag model [s^ as our approximation to baryon structure, we note that there has 
been a parallel development [s^ based upon the covariant, chiral symmetric NJL model [33 |. with quark confinement 
modelled using the proper time regularization proposed by the Tiibingen group [H, [s^ . The latter model has many 
advantages for the computation of the medium modification of form factors and structure functions, with the results 
for spin structure functions [13, HI] offering a unique opportunity to test the fundamental idea of the QMC model 
experimentally. However, in both models it is the effect of quark confinement that leads to a positive polarizability 
and a natural saturation mechanism. 

Although the underlying physics of QHD and QMC is rather different, at the hadronic level the equations to be 
solved are very similar. We therefore focus on the changes which are required. 

1. Because of the scalar polarizability of the hadrons, which accounts for the self-consistent response of the 
internal quark structure of the baryon to the applied scalar field jl^, the effective masses appearing in QMC are 
non-linear in the mean a field. We write them in the general form 



Ms = Mb - w% 9aN{o) + i^'^B {9<yN{(^) f , 



(24) 
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FIG. 1: (Color online) Baryon effective masses within the QMC model, parameterized as a function of the mean scalar field, 
(a). The values at (a) = are the vacuum masses as found in Table|T] We show the effective masses only up to (cr) = 100 MeV 
which corresponds to about 2 fm"'' (6-8 po), beyond which higher order terms not shown in Eq. (|24p become significant. 



where the weightings, w^, Wg, and the scalar polarizability of the nucleon, rf, must be calculated from the underlying 
quark model. Note now that only the coupling to the nucleons, gaN, is required to determine all the effective masses. 

The most recent calculation of these effective masses, including the in-medium dependence of the spin dependent 
hyperfine interaction yields the explicit expressions: 



Mjv((a)) 
Ma ((a)) 



0.0022 + 0.1055i?5^°° - 0.0178 {R%'"'f] {gaN{<j)f , 



Ma 



0.6672 + 0M62R%'"' ~ 0.0021 {R%^'') g^N (ff) 



0.0016 + 0.0686i?^S=° - 0.0084 (R%'"' 



{9<yN{<y)) , 



Me((ct)) = Ms- 



0.6706 - 0.0638i? 



0.008 (R 



frcc\ 



9aN{<y) 



(25) 



-0.0007 + 0.0786i?^'''= - 0.0181 {R%''''f (.9<Tiv(cr))^ , 



MH((a)) 



Me 



0.3395 + 0.02822i?5;/°'= - 0.0128 {R 



9(yN (cr) 



-0.0014 + 0.0416i?^'=<= - 0.0061 (i?^°°) (5<xw(o-))' 



We take R^^'^ = 0.8 fm as the preferred value of the free nucleon radius, although in practice the numerical results 
depend only very weakly on this parameter p^ . 

Given the parameters in Eq. I|25p . all the effective masses for the baryon octet are entirely determined. They are 
plotted as functions of (cr) in Fig. [1] and we see clearly that they never become negative. (Note that the range of (cr) 
covered here corresponds to densities up to (6-8)po)- 

2. Since the mean scalar field, (cr), is derived self-consistently by taking the derivative of the energy density with 
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respect to (a), the scalar field equation 



^ ' ^ ml (27r)3 7 ^fc2 + (Af* )2 ^ ^ 



has an extra factor, denoted by 



C((<7)) = K-^Sdg,^(a)]. (27) 



Note that the d term (the scalar polarizability) in C{{a)) does not have the factor of i that is found Eq. (|24p . because 
of the differentiation. 

Given this new term in the equation for the mean scalar field, we can see that this allows feedback of the scalar 
field which is modelling the internal degrees of freedom of the baryons. This feedback prevents certain values of (cr) 
from being accessed. 

3. The couplings to the proton are re-determined by the fit to saturation properties (minimum binding energy 
per baryon and saturation density) with the new effective masses for the proton and neutron. The couplings for QMC 
which provide a fit to saturated nuclear matter are shown in Table IIIII 

Given these changes alone, QHD is transformed into QMC. When we compare the results of Section fVlII with those 
of Ref. [l^ minor differences arise because the QMC calculations in Ref. [3 are performed at Hartree-Fock level, 
whereas here they have been performed at Hartree level (mean-field) only. 

TABLE III: Couplings for QMC with the octet of baryons, fit to saturation of nuclear matter. 







9p 


8.278 


8.417 


8.333 



IV. DECONFINED QUARK MATTER 

We consider two models for a deconfined quark matter phase, both of which model free quarks in /3-equilibrium. The 
first model, the MIT bag model [121, commonly used to describe the quark matter phase because of its simplicity. 

In this model we consider three quarks with fixed masses to possess chemical potentials related to the independent 
chemical potentials of Eq. pT|) via 

12 11 

= - gMe, A'd = + gA^e, = Md, (28) 

where quarks have a baryon charge of i since baryons contain 3 quarks. Because the quarks are taken to be free, the 
chemical potential has no vector interaction terms, and thus 



Hq =^ J Pp + m'i ; q€{u,d,s}. (29) 



The EoS can therefore be solved under the conditions of Eq. ([23]) . 

As an alternative model for deconfined quark matter, we consider a simplified Nambu-Jona-Lasinio (NJL) model [s^ . 
in which the quarks have dynamically generated masses, ranging from constituent quark masses at low densities to 
current quark masses at high densities. The equation for a quark condensate at a given density (and hence, kp) in 
NJL is similar to the scalar field in QHD/QMC, and is written as 

f 1 A/; eikp - \k\) 9(A - kp) , 

J (2^^ ) Jk-^ + {M*y 

where M* denotes the kp dependent (hence, density dependent) quark mass; A/'c is the number of color degrees of 
freedom of quarks; and A is the momentum cutoff. This is self-consistently calculated via 



M* =Mcurrent-G(V'V'), 



(31) 
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where G is the couphng and Mcurront the current quark mass. 

To solve for the quark mass at each density, we must first find the couphng, G, which yields the required constituent 
quark mass in free space {kp — 0). The coupling is assumed to remain constant as the density rises. In free space, 
we can solve the above equations to find the coupling 



G = 



[m; - A/c, 



1 M*„ 



\k\-kF)e{K-kF) 3 



(27r)3 



fc2 



(Tk 



(32) 



We solve Eqs. ([5DH5^ for Nc = 3 to obtain constituent quark masses of Mu,d = 350 MeV using current quark masses 
of Mcurrent = 10 McV for the light quarks, and to obtain a constituent quark mass of Mg = 450 MeV using a current 
quark mass of Mcurrent = 160 McV for the strange quark, with a momentum cutoff of A = 1 GcV. Ai kp ~ wc find 
the couplings to be 



G^,d = 0.148 fm^ G, = 0.105 fm^ 



(33) 



We can now use these parameters to evaluate the dynamic quark mass M*, for varying values of kp, by solving 
Eq. (j30p and Eq. ([3T|) self-consistently. The resulting density dependence of M* is illustrated in Fig. [2l This shows 
that the masses of the quarks eventually saturate and are somewhat constant above a certain density. We can then 
construct the EoS in the same way as wc did for the MIT bag model, but with density-dependent masses, rather than 
fixed masses. 



500 




1 2 3 4 5 6 

[fm-^] 

FIG. 2: (Color online) Density dependent (dynamic) masses for quarks using NJL. The mass &t kp = Q is the constituent 
quark mass, and the mass at the cutoff oi kp = K is roughly the current quark mass. This model successfully reproduces the 
behaviour found within the Schwinger-Dyson formalism, and we consider the model to be more sophisticated than the constant 
quark mass MIT bag model. 
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V. PHASE TRANSITIONS 
A. Equilibrium Conditions 



We now have a description of hadronic matter with quark degrees of freedom, but we are still faced with the issue 
that the baryons are very densely packed. We wish to know if it is more energetically favourable for deconfined quark 
matter to be the dominant phase at a certain density. To do this, we need to find a point (if it exists) at which 
stability is achieved between the hadronic phase and the quark phase. 

The condition for stability is that chemical, thermal, and mechanical equilibrium between the hadronic (H) and 
quark (Q) phases is achieved, and thus that the independent quantities in each phase are separately equal. Thus 
the two independent chemical potentials, {fin,fJ-e), are each separately equal to their counterparts in the other phase, 
i.e. {^in)H ~ {^J'n)Q, and {ne)H = (Me)Q (chcmical equilibrium); the temperatures are equal {Th ~ Tq) (thermal 
equilibrium); and the pressures are equal {Ph = Pq) (mechanical equilibrium). For a discussion of this condition, see 
Ref. [3^. We consider both phases to be cold on the nuclear scale, and assume T = 0, so the temperatures are by 
construction equal. We must therefore find the point at which, for a given pair of independent chemical potentials, 
the pressures in both the hadronic phase and the quark phase arc the same. 

To find the partial pressure of any baryon, quark, or lepton species, i, we use 

3(27r)3 J ^k^ + (M*)2 ^ ^ 

where TVc = 3 for quarks, and TVc = 1 for baryons and leptons. To find the total pressure in each phase we use 

Ph = Y.Pb+Y.P^+Y. \<{-^f-\<{<^?^ (35) 

B t OL—LJ^p 

which is equivalent to Eq. (fT7|) . and 

q I 

where B in the quark pressure is the bag energy density. For the QMC model described in Section HTll and a Fermi gas 
of quarks, both with interactions with leptons for charge neutrality, a point exists at which the condition of stability, 
as described above, is satisfied. 

At this point, it is equally favourable that hadronic matter and quark matter are the dominant phase. Beyond this 
point, the quark pressure is greater than the hadronic pressure, and so the quark phase has a lower thermodynamic 
potential (through the relation P = —£7) and the quark phase will be more energetically favourable. To determine 
the EoS beyond this point, we need to consider a mixed phase. 



B. Mixed Phase 



We can model a mixed phase of hadronic and quark matter — as opposed to modelling a simple direct phase 
transition between the two, a Maxwell construction, which would have a discontinuity in the density, while retaining a 
constant pressure between the two phases — using the method of Glendenning. A detailed description of this appears 
in Ref. [11]. 

We solve for the hadronic EoS using the independent chemical potentials as inputs for the quark matter EoS, as the 
order parameter, p, the conserved baryon density, increases until we find a point (if it exists) at which the pressure 
in the quark phase is equal to that of the hadronic phase. Once we have the density and pressure at which the phase 
transition occurs, we change the order parameter from the conserved baryon density to the quark fraction, x- If 
consider the mixed phase to be a fraction of the hadronic matter and a fraction of the quark matter, then the mixed 
phase (MP) of matter will have the following properties; the total density will be 

PMP = (1 - X) PHP + X PQP, (37) 

where php and pqp are the densities in the hadronic and quark phases, respectively. The equivalent baryon density 
in the quark phase. 



^QP = X! ^9 " 3(p« + + Ps)-, 
q 



(38) 
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arises because of the restriction that a bag must contain 3 quarks. 

According to the condition of mechanical equihbrium, the pressure in the mixed phase wiU be 

Pmp = Pup = Pqp- (39) 

We can step through values < x < 1 and find the density at which equilibrium is achieved, keeping the mechanical 
stability conditions as they were above. In the mixed phase we need to alter our definition of charge neutrality; it 
becomes possible now that one phase is (locally) charged, while the other phase carries the opposite charge, making 
the system globally charge neutral. This is achieved by enforcing 

= {1-X) PHP + X Pqp + Pi, (40) 

where this time we are considering charge densities, which are simply charge proportions of density and is the 
lepton charge density. For example, the charge density in the quark phase is given by 

^ 2 11 

We continue to calculate the densities until we reach x = 1, at which point the mixed phase is now entirely charge 
neutral quark matter. After this point, we continue with the EoS for pure charge neutral quark matter, using p as 
the order parameter. 



Pressure [MeV fm ] 

Quark 
Hadron 
Projection 

600 n 




FIG. 3: (Color online) Illustrative locus of values for /ie,/in,-P for phases of hadronic matter and deconfined quark matter. 
Note that pressure increases with density, and that a projection onto the UnfJ-e plane is a single line, as ensured by the chemical 
equilibrium condition. 
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VI. STELLAR SOLUTIONS 



To test the predictions of these models, we find solutions of the Tolman-Oppenheimer-Volkoff (TOV) [40| equation 

dP _ G{P + £) {M{R) + AttR^P) 



dR R{R-2GM{R)) 
where the mass, M{R), contained within a radius R is found by integrating the energy density 



(42) 



M{R) = i Anr^S dr, (43) 

and £ and P are the energy density and pressure in the EoS, respectively. 

Given an EoS and a choice for the central density of the star, this provides static, spherically symmetric, non- 
rotating, gravitationally stable stellar solutions for the total mass and radius of a star. For studies of the effect of 
rapid rotation in General Relativity we refer to Refs. [4ll.[43|. This becomes important for comparison to experimental 
data, as only data for stellar masses exists (with the single, disputed exception from [T^l. we can use the model to 
predict the radii of the observed stars. 



VII. RESULTS 



To obtain numerical results, we solve the meson field equations, Eqs. ([7]-[9]), with the conditions of charge neutrality, 
fixed baryon density, and the equivalence of chemical potentials given by Eq. (|23p . for various models. Having found 
the EoS by evaluating the energy density, Eq. (fTBl) . and pressure, Eq. ([T7]), we can solve for stellar solutions for an 
EoS using the TOV equation. The radius of the star is defined as the radius at which the pressure is zero and is 
calculated using a fourth-order Runge-Kutta integration method. 

The EoS for octet QMC hadronic matter is shown in Fig. 2] alongside the same model when including a phase 
transition to 3-flavor quark matter modelled with the MIT bag model, and the results do not appear to differ much at 
this scale. The theoretical causality limit of P = £ is also shown (corresponding to the limit Wsound = c) and we can 
see that these models do not approach this limit at the scale displayed. This is because of the softening of the EoS 
that occurs with the introduction of hypcrons, enlarging the Fermi sea to be filled and reducing the overall pressure. 

The species fraction for each particle, Yi, is simply the density fraction of that particle, and is calculated via 

= iG{p,n,A,I]+,I]",I]-,S°,S-,^,g}, (44) 

P 

where p is the total baryon density. The species fractions for octet QMC when a phase transition is neglected are 
shown in Fig. [5l where we note that the A species fraction is enhanced and the E species fractions are suppressed 
with increasing density. From the investigations by Rikovska-Stonc et al. fl6l | we expect that the S would disappear 
entirely if we were to include Fock terms. 

The value of compression modulus and effective nucleon mass at saturation are frequently used as a comparison to 
experimental evidence. Models which neglect quark-level interactions, such as QHD, typically predict much higher 
values for the compression modulus than experiments suggest. In the symmetric (nuclear) matter QHD model 
described in this paper, we find values of (M*/M)sat = 0.56 and K = 525 MeV which are in agreement with p^ . 
but as stated in that reference, not with experiment. For QMC we find a significant improvement in the compression 
modulus; K = 280 MeV which lies at the upper end of the experimental range. The nucleon effective mass at 
saturation for QMC is found to be (M*)sat = 735 MeV, producing (M*/M)sat = 0.78. 

When we calculate the EoS including a mixed phase and subsequent pure quark phase, we find that small changes 
in the parameters can sometimes lead to very significant changes. In particular, the bag energy density, B, and 
the quark masses in the MIT bag model have the ability to both move the phase transition points, and to vary the 
constituents of the mixed phase. Wc have investigated the range of parameters which yield a transition to a mixed 
phase and these are summarised in Table IIVI For illustrative purposes we show an example of species fractions for 
a reasonable set of parameters [B'^/^ = 180 MeV and mu,d,s = 3,7,95 MeV) in Fig. [51 Note that in this case the A 
hyperon enters the mixed phase briefly (and at a low species fraction). 

Note that the transition density of pmp ~ 0.12 fm~'^ produced by the combination of the octet QMC and MIT bag 
models (as shown in Fig. [6]) is clearly not physical as it implies the presence of deconfined quarks at densities less 
than po- 

With small changes to parameters, such as those used to produce Fig. [7] in which the bag energy density is given 
a slightly higher value from that used in Fig. [S] (S^/^ increased from 180 MeV to 195 MeV, but the quark masses 
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FIG. 4: (Color online) Equation of State for; nucleonic 'N' matter modelled with octet QMC but where hyperons are explicitly 
forbidden; nucleonic matter where a phase transition to NJL modelled quark matter is permitted; baryonic 'N+Y' matter 
modelled with octet QMC including hyperons; and baryonic matter where a phase transition to MIT bag modelled quark 
matter is permitted. The line P = £ represents the causal limit, iigound = c. The bends in these curves indicate a change in the 
composition of the EoS, such as the creation of hyperons or a transition to a mixed or quark phase. Note that at low energies 
(densities) the curves are identical, where only nucleonic matter in /3-equilibrium is present. 



remain the same), it becomes possible for the 5 hyperons to also enter the mixed phase, albeit in that case with small 
species fractions, Y^,Ys < 0.02. 

The TOV solutions for octet QMC with and without a phase transition to a mixed phase are shown in Fig. [S] 
The stellar masses produced using these methods are similar to observed neutron star masses. Once we have solved 
the TOV equations, we can examine individual solutions and determine the species content for specific stars. If we 
examine the solutions with a stellar mass of M = 1.2 Mq, where Mq is a solar mass, for the set of parameters used to 
produce Figs. [5] and [51 wc can find the species fraction as a function of stellar radius to obtain a cross-section of the 
star. This is shown in Fig. [9] for the case of no phase transition, and Fig. [10] for the case where we allow a transition 
to a mixed phase, and subsequently to a quark matter phase. 

If we now examine the stellar solution with mass M = 1.2 A/q of the set of parameters (_B^/^ = 195 MeV and 
mu,d,s = 3, 7, 95 MeV used to produce Fig. [7]) as shown in Fig. [11] we note that the quark content of this 10.5 km star 
reaches out to around 8 km, and that the core of the star contains roughly equal proportions of protons, neutrons 
and A hyperons with Yi ~ 10%. 

Within a mixed phase, we require that for a given pair of fin and fie at any value of the mixing parameter x, the 
quark density is greater than the hadronic density. This condition ensures that the total baryon density increases 
monotonically within the range pQp > pMP > PHP, as can be seen in Eq. (|37p . An example of this is illustrated in 
Fig. [m for a mixed phase of octet QMC and 3-flavor quark matter modelled with the MIT bag model. 

The use of quark masses corresponding to the NJL model results in a quark density that is lower than the hadronic 
density, and as a result there are no solutions for a mixed phase in which the proportion of quarks increases with 
fraction x, while at the same time the total baryon density increases. It may be possible that with smaller constituent 
quark masses at low density, the Fermi momenta would provide sufficiently high quark densities, but wc feci that it 
would be unphysical to use any smaller constituent quark masses. This result implies that, at least for the model we 
have investigated, dynamical chiral symmetry breaking (in the production of constituent quark masses at low density) 
prevents a phase transition from a hadronic phase to a mixed phase involving quarks. 
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FIG. 5: Species fractions, Yi, for octet QMC where a transition to a mixed phase is exphcitly forbidden. Note that in this 
case, all of the octet baryons contribute at some density, and that with increasing density the species fractions of E hyperons 
are suppressed while the A species fraction is enhanced. Parameters used here are shown in Table Hill 



We do note, however, that if we restrict consideration to nucleons only within the QMC model (with the same 
parameters as octet QMC), and represent quark matter with the NJL model, we do in fact find a possible mixed 
phase. More surprisingly, the phase transition density for this combination is significantly larger than the case where 
hyperons are present. An example of this is shown in Fig. 1131 with parameters found in Table ITVl This produces 
a mixed phase at about 3po {p = 0.47 fm""^) and a pure quark matter phase above about 10.5po (p = 1-67 fm~'^). 
We note the coincidence of this phase transition density with the density corresponding to one nucleon per nuclcon 
volume, with the aforementioned assumption of R^^'^ = 0.8 fm, though we do not draw any conclusions from this. 
Performing this calculation with quark matter modelled with the MIT bag model produces results similar to those of 
Fig. [6] except of course lacking the A hyperon contribution. Although this example does show a phase transition, the 
omission of hyperons is certainly unrealistic. This docs however illustrate the importance and significance of including 
hyperons, in that their inclusion alters the chemical potentials which satisfy the equilibrium conditions in such a way 
that the mixed phase is no longer produced. 

For each of the cases where we find a phase transition from baryonic matter to quark matter, the solution consists 
of negatively charged quark matter, positively charged hadronic matter, and a small proportion of leptons, to produce 
globally charge neutral matter. The proportions of hadronic, leptonic and quark matter throughout the mixed phase 
(for example, during a transition from octet QMC matter to 3-flavor quark matter modelled with the MIT bag model) 
are displayed in Fig. [Ml A summary of the results of interest is given in Table IIVI 

Results for larger quark masses are not shown, as they require a much lower bag energy density to satisfy the 
equilibrium conditions. For constituent quark masses, we find that no phase transition is possible for any value of 
the bag energy density, as the quark pressure does not rise sufficiently fast to overcome the hadronic pressure. This 
is merely because the mass of the quarks does not allow a sufficiently large Fermi momentum at a given chemical 
potential, according to Eq. ([29)) . 



VIII. CONCLUSIONS 



We have produced several EoS that simulate a phase transition from octet QMC modelled hadronic matter, via 
a continuous Glendenning style mixed phase to a pure, deconfined quark matter phase. This should correspond to 
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FIG. 6: Species fractions, Yi, for octet QMC (the same as in Fig. [5]) but where now we allow the phase transition to a mixed 
phase involving quark matter modelled with the MIT bag model, and subsequently to a pure deconfined quark matter phase. 
Parameters used here are summarised in Table IIVI Note that with these parameters, the A is the only hyperon to appear in 
the mixed phase, and does so at a much higher density than the case where the transition to a mixed phase is forbidden. We 
also note that with these parameters, the transition to a mixed phase occurs below saturation density, po- 



a reasonable description of the relevant degrees of freedom in each density region. The models used here for quark 
matter provide a framework for exploring the way that this form of matter may behave, in particular under extreme 
conditions. The success of the QMC model in reproducing a broad range of experimental data gives us considerable 
confidence in this aspect of these calculations, and provides a reasonable hadronic sector and calculation framework, 
which then awaits improvement in the quark sector to produce realistic stellar solutions. 

We have presented EoS and stellar solutions for octet QMC matter at Hartree level. We have explored several 
possible phase transitions from this hadronic sector to a mixed phase involving 3-flavor quark matter. The corre- 
sponding EoS demonstrate the complexity and intricacy of the solutions as well as the dependence on small changes 



TABLE IV: Table of species content {N = nucleons, Y = hyperons, I = leptons, q = quarks); inputs {B^^^, mq); and results 
for octet QMC and quark models presented in this paper. pY, Pmp and pqp represent the density at which hyperons first 
appear (A is invariably the first to enter in these calculations); the density at which the mixed phase begins; and the density 
at which the quark phase begins, respectively. Figures for selected parameter sets are referenced in the final column. Dynamic 
NJL quark masses are determined by Eqs. (|30ft32|l . 

Particles: B^/* (MeV) {mu,md,ms} (MeV) py (fm"'^) pMP (fm^^*) pqp (fm^^) Figure: 

N, Y, ^ — — 0.27 — — Fig. [5] 

N, Y, ^, q 180 {3,7,95} 0.55 0.12 0.95 Fig. |6] 

N, Y, ^, q 195 {3,7,95} 0.35 0.24 1.46 Fig. [7] 

N, Y, £, q 170 {30, 70, 150} 0.56 0.10 0.87 — 

N, Y, £, q 175 {100, 100, 150} 0.44 0.16 1.41 — 

N, ^, q 180 Dynamic (NJL) — 0.47 1.67 Fig. [131 
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FIG. 7: Species fractions, Yi, for octet QMC (the same as in Fig. [6] but now where the bag energy density has been increased 
to B^^^ = 195 MeV). Note that now tlie appearance of hyperons occurs at a smaller density than in the case of Fig. [6l the 
transition to a mixed phase occurs at a slightly larger density, and that now H hyperons are present in the mixed phase. 



in parameters. The stellar solutions provide overlap with the lower end of the experimentally acceptable range. 

Several investigations were made of the response of the model to a more sophisticated treatment of the quark masses 
in-medium, namely the NJL model. In that model the quark masses arise from dynamical chiral symmetry breaking 
and thus take values typical of constituent quarks at low density and drop to current quark masses at higher densities. 
The result is that no transition to a mixed phase is possible in this case. 

The omission of hyperons in the QMC model yields a transition to a mixed phase of cither NJL or MIT bag model 
quark matter, as the hadronic EoS is no longer as soft. This observation makes it clear that hyperons can play a 
significant role in the EoS. However, we acknowledge that their presence in neutron stars remains speculative. 

The models considered here reveal some important things about the possible nature of the dense nuclear matter 
in a neutron star. It seems that if dynamical chiral symmetry does indeed result in typical constituent quark masses 
in low density quark matter, then a phase transition from hadronic matter to quark matter is unlikely. This result 
invites further investigation. 

The results presented in Fig. [8] indicate that the model in its current form is unable to reproduce sufficiently massive 
neutron stars to account for all observations, notably the observed stellar masses of 1.45 Mq and larger. This is a 
direct result of the softness of the EoS. This issue will be explored in a future publication via the inclusion of Fock 
terms, which have been shown to have an effect on the scalar and vector potentials 43 1. 

Many open questions remain to be investigated in further work, including the effects of Fock terms, and the density 
dependence of the bag energy density in the quark phase, which can be calculated explicitly within the NJL model. 
The quark matter models used here are still not the most sophisticated models available, and further work may involve 
an investigation of the effects of color-superconducting quark matter 0, [4^ . 
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FIG. 9: Species fractions for octet QMC in /3-equilibrium, where the phase transition to a mixed phase is explicitly forbidden, 
as a function of stellar radius for a stellar solution with a total mass of 1.2 Af©. The parameters used here are the same as 
those used to produce Fig. [S] 
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FIG. 10: Species fractions for octet QMC with a phase transition to 3-flavor quark matter modelled with the MIT bag model, 
as a function of stellar radius for a stellar solution with a total mass of 1.2 Mq. The parameters used here are the same as 
those used to produce Fig. |6l Note that in this case one finds pure deconfined 3-fiavor quark matter at the core (all of some 
3.5 km) of this star, and still a small proportion of A in the mixed phase. 
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FIG. 11: Example of the interior of a star of total stellar mass M — 1.2Mq where the bag energy density is given a slightly 
higher value from that used in Fig. 1101 (increased from B^''* — 180 MeV to 195 MeV), but the quark masses remain the same. 
This illustrates that with relatively minor adjustments to the parameters, large changes can be introduced to the final solution. 
In this case H hyperons can provide a nonzero contribution to the composition of a star. Note that in this case, quark matter 
appears at 8 km, and at the core there exists a mixed phase containing nucleons, quarks, as well as A, EP, and H~ hyperons. 
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FIG. 12: (Color online) Densities in the mixed phase for octet QMC mixed with 3-flavor quark matter modelled with the MIT 
bag model. Note that at all values of x (the mixing parameter according to Eq. (jSTjl), the equivalent quark baryon density 
is greater than the hadronic baryon density, allowing the total baryon density to increase monotonically. The total density is 
found via Eq. ((37)). 
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FIG. 13: Species fractions for a phase transition from QMC nuclear matter to 3-flavor quark matter modelled with NJL. 
Note that in this unphysical case, a phase transition is possible, and occurs at a value of about p — 0.47 fm~"^. We note the 
coincidence with the density of one baryon per baryon volume, but draw no conclusions from this. A similar transition from 
QMC nuclear matter to 3-flavor quark matter modelled with the MIT bag model produces results almost identical to those of 
Fig. [6] except that in that case there is no contribution from the A hyperon. 
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FIG. 14: (Color online) Charge densities (in units of the proton charge per cubic fm) in the mixed phase for a transition from 
octet QMC to 3-fiavor quark matter modelled with the MIT bag model. Note that following the mixed phase, the quarks are 
able to satisfy charge neutrality with no leptons. x is the mixing parameter within the mixed phase according to Eq. p7|l . 



